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Starting from a set of coupled Boltzmann equations, we investigate the thermalization of a two- 
species cold atomic gas confined either in a box or in an isotropic harmonic trap. We show that the 
thermalization times, by contrast to the collision rate, depend on the interferences between scattering 
partial waves. The dynamics of thermalization in a harmonic trap is also strongly dependent upon 
the ratio between the collision rate and the trap frequencies. 
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I. INTRODUCTION 



Ultracold mixtures of alkali atoms have attracted great 
interest among atomic physicists. The ability to simul- 
taneously cool and trap multiple species offers the pos- 
sibility to study a large variety of physical phenomena 
not previously accessible with a single species such as 
novel degenerate quantum systems made of Fermi-Bosc 
mixtures 

Thermalization plays a crucial role in the cooling of 
a sample by evaporation, and/or in the implementation 
of sympathetic cooling. The detailed knowledge of the 
underlying physics is all the more important for exper- 
iments where the time available to perform evaporation 
is limited [2J , or for the optimization of the evaporation 
ramp in the multi-partial- wave regime [3]. In addition 
the return to equilibrium of a gas initially prepared in an 
out-of-equilibrium state is a widely used way to infer the 
value of the scattering length 0, [E H, 0, ll.l9l.IToj. 

The purpose of this paper is to provide a theoretical 
framework to describe the thermalization between two 
cold gases confined in a box or a harmonic trap. In 
contrast to previous theoretical studies on this subject 
@, S EH ! we explicitly take into account all partial waves 
involved in the collision process, and we work out the role 
played by the confinement in all collisional regimes. Since 
the basic mechanism described by the Boltzmann equa- 
tion is binary collisions, two-component systems exhibit 
all the complexities of n-component systems. The gener- 
alization of the results we obtain to systems containing 
more than two species is therefore straightforward. 

The paper is organized as follows. Section [TT1 describes 
the theoretical framework used to investigate the dynam- 
ics of thermalization. Section Hill is devoted to the simple 
case of a gas confined in a box. It gives the explicit ex- 
pression for the rethermalization time including the in- 
terference terms between all partial waves involved in the 
collision process. Section HVl addresses the same problem 
for a gas trapped in an isotropic harmonic trap. 



II. THEORETICAL FRAMEWORK 

The gas i (i — 1,2) consists of 2Vj particles of mass rrii 
initially thermalized at Tf . The dynamics of the gas is 
described by a set of two coupled Boltzmann equations 
for each phase-space distribution function /j(r, Vj,i): 

df 

-g + ifaHi} = Iu[fi]+Iij[fiJj], i¥>J- (1) 

The distribution function fi for the species i is gov- 
erned by the hamiltonian Hi through the Poisson bracket 
{fi, Hi}, and by binary collisions. The appearance of two 
collision terms on the right-hand side of Eq. dTJ is a re- 
sult of the possibility of fa to change through both self 
collisions (i-i collisions) or cross collisions (1-2). The ex- 
pression for the collision integrals is 

WiJj] = I d 2 n>d 3 v B \v B -v A \^ x 

[fi (ry A , t)fj (ry B) t)- fi (r, v A) t)fj (r, v B , *)], (2) 

where i,j = 1,2. They account for elastic collisions be- 
tween particles labelled A and B. We denote the ve- 
locities of the two atoms before they interact with each 
other by and yg, and after the collision by v'^ and 
v'b- The expression for the collision terms can be easily 
extended to include the effects of both Bose and Fermi 
statistics leading to the Boltzmann-Nordheim equation 
El- In fact, most of the results discussed in this 
paper also hold in the presence of quantum degeneracy, 
provided the system is not Bose-Einstein condensed. The 
quantitative estimates of collision effects presented in this 
work are, however, based on classical statistics. As shown 
in Refs [ll| for bosons and [l4| for fermions, our ap- 
proach, based on the result of the classical physics, is 
valid as soon as the temperature of the sample is larger 
than twice the temperature for degeneracy, i.e. the crit- 
ical temperature for bosons and the Fermi temperature 
for fermions. 

The collision problem is simplified by an appropriate 
change of variables. We introduce the center of mass ve- 
locity vq = (mAVA + TriB'VB)/('mA + rnB), the total mass 
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M = rriA + tub, the relative velocity v r = - vb, 
and the reduced mass /i = m^m^/M. The relative ve- 
locity changes from the value v r before the collision to 
the value v' r after the collision, keeping the same value 
for the modulus v r = v' r because of energy conservation. 
The collision process is described by specifying the scat- 
tering direction with the polar angle 9' and the azimuthal 
angle ip 1 of the final relative velocity v' r with respect to 
the initial relative velocity v r before the collision. The 
corresponding solid angle d 2 Q! = d<p'sm9'd9' enters the 
expression for the integrand of the collision integrals. 

The angular dependence of the elastic scattering di- 
rectly derives from the partial-waves decomposition of 
the scattering amplitude in the quantum mechanical de- 
scription of interactions (l5j . For non-identical atoms, 
the asymptotic form of the scattering amplitude f(k, 0') 
reads 

2=0 

where k — /j,v r /Ti, 0' <E [0, 7r) and P;[cos(0')] are the Leg- 
endre polynomials. All the complexity of the interatomic 
potential is contained in the fc-dependence of the phase 
shifts 5i. The differential cross section for non-identical 
atoms is given by da^/dQ, 1 = |/(fc,6>')| 2 . 

For identical atoms, we have to take into account the 
(anti)symmetrization principle. Accordingly the differ- 
ential cross section takes the form: 

^ = \ f (k,9') + ef(k,n-9')\ 2 , (4) 

where e = f for bosons and e = — 1 for fcrmions. The 
integration must be carried out only in the half sphere 
9' e [0, tt/2). As a consequence of the parity of the spher- 
ical harmonic, the only partial waves contributing to 
the scattering cross-section for polarized bosons (respec- 
tively, fermions) correspond to even (respectively, odd) 
values of I. The interference between partial waves are 
contained in the differential cross section [l6l |. 

A relevant parameter to describe the thermalization of 
the gases is the collision rate. The average total number 
of collisions per unit of time for atom species i and 
j is obtained by integrating over space and velocity the 
output channel term of the collision integral in Eq. : 

Tij = J d 3 rd 3 VAd 3 v B (Tij\vA-VB\fi(r,VA,t)fj(r,VB,t) 

From this quantity we deduce the expression for the mean 
collision rate per atom 7^ = Tij/Ni. The integrated 
cross section Uij is just the sum of the contributions from 
partial-waves: 

(Tij K) = (1 + Sij) V —{21 + 1) sin 2 Si, (6) 
1 k 

where the factor (1 + <5j j) accounts for the constructive 
interference of undistinguishable scattering channels for 



identical atoms, 8ij being the Kronecker delta symbol. 
The sum in Eq. ^ is taken over even (resp. odd) values 
for identical (i = j) bosons (resp. fermions), and over 
all integer values for non-identical atoms (i j). The 
expression for the collision rate is just a simple sum over 
the partial waves, it does not exhibit interference. This is 
to be contrasted with the expression of the thermalization 
rate as will be illustrated in the next section. 

The determination of the thermalization rate is based 
on the method of averages combined with an appro- 
priate ansatz. Starting from Eq. |T]), one readily derives 
the equation for the average value of a general dynamical 
quantity C(r, v*): 

^ + (0{fi, = {OIn)i + (Ohiii, (7) 

where the average is taken in both position and velocity 
space: (0)i — J d 3 rd 3 Vifi(r,v i7 t)0(r,Vi)/Ni. As a con- 
sequence of the collisional invariants — number of atoms, 
momentum, and energy — , {OiIa)i = if Oi is of the 
form L = a(r) + b(r) ■ v, + c(r)v 2 . 



III. THERMALIZATION IN A BOX 

As a first application of the model, we consider the 
thermalization of a two-species gas in a box of volume 
V. We assume that atoms can undergo only specular 
reflection on the walls of the box, which can be realized 
experimentally by using the dipolar optical force [l8l |. 
Initially, the gases are at different temperatures T-f 7^ T$. 
After thermalization, the temperature of the gases will 
be the same T f = (Nxlf + N 2 T°)/(N 1 + N 2 ), within 
the Boltzmann dynamics. To evaluate the thermalization 
time, we write the equation for the mean total kinetic 
energy (A^miV 2 i/2)j of species i: 

The Poisson bracket of Eq. ([7]) vanishes for a homoge- 
neous system, and the self collision integral term does not 
contribute because of conservation of kinetic energy for 
elastic collisions. The calculation of the collision integral 
term requires an ansatz for the distribution function. We 
choose a Gaussian ansatz for the phase-space distribution 
function of each species: 

where M% = Nim 3 ' 2 /(V(2nkBTi) 3 ^ 2 ) is the normaliza- 
tion factor. The time dependence is contained only in 
the effective temperatures Tj. We assume that (T° — 
T f )/T f < 1 and (T 2 ° - T } )/T } < 1. Consequently, we 
can evaluate the total number of cross collisions per unit 
of time by setting T x = T 2 = T f in Eq. © with the 
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Gaussian ansatz for the distribution functions: 

Tl2 = —j=— TJ—C (T12), 
\/7T V 



(10) 



where we have introduced the velocity c = (2ksTf / ix) 1 / 2 
and the thermally averaged cross section {(J12) — 
2 / °° dxon{cx)x z e~ x . 

The expansion around the final temperature Tf of the 
collision terms yields a set of linear equations. We can 
rewrite them in the form (see Appendix IX)) : 



d{T, - Tj 



dt 



Ni + N< \ Ti - T, 



with 



8M «^ia) 
3M (a-ia) 



12, 



(11) 



(12) 



where (((T12)) = dx ai2{cx)x 5 e~ 
d\2(cx) is defined as: 

012 (ex) = 2vr / sin 0' (1 - cos 0' 
Jo 



The quantity 



dfl' 



(13) 



The new feature of this expression arises from the inter- 
ference between partial waves resulting from the angular 
integration. This is in contrast to the expression for the 
total cross section, which just contains a simple sum over 
partial waves. We show in Appendix [Al that interference 
occur only between partial waves differing by at most one 
unit of angular momentum. 

We thus find that the relaxation corresponds to an 
exponential decay but the thermalization times r is not 
in general proportional to the collision rate because it 
contains partial wave interferences terms. 

At very low temperature where only s- waves contribute 
to the collision process, no interference can occur and the 
relaxation time is proportional to the inverse of the col- 
lision rate. If the cross section depends on the relative 
wave vector k of the collision, the factor of proportional- 
ity depends on temperature. 

For a cross section of the form <7i2(fc) = 47ra^ 2 /(l + 
a 2 2 k 2 ), one finds 



(012) 



4™? 2 £[1 
27m? 2 £[l 



Vr(o,o], 
£ + £Vr(o,0], 



(14) 



where £ = h 2 /(2fik B T f a 2 12 ) and T(a,z) = J™ t^e^dt. 
In the very low temperature limit £ 3> 1, we have 
(012) = ((< T i2)) = 47ra^ 2 ; in the unitary limit ( < 1, 
we obtain (012) = 2((<7i 2 )) = 2ttT% 2 / {^ksTf). As a con- 
sequence, the number (7T12) of inter-species collisions 
per atom required to equilibrate the temperature of a 
two-component system made of atoms of the same mass 
varies from 1.5 in the very low temperature limit to 3 in 
the unitary limit according to Eq. (TT2")) . 

In the following, we explicitly derive the curves for the 
thermalization time r and the inter-species collision rate 



7 for the same number of 87 Rb atoms N\ — N2 = N 
in two different internal states \5Si/ 2 ,F = 2,m = 1) 
and 155*1/2, F — l,m = —1). This calculation has been 
performed by taking into account the first five partial 
wave contributions for the cross collisions [lj|. We plot 
in Fig. (p} the average cross section ((a 12)) normalized to 
its value with the cross section taken at zero energy. The 
resonance shape results from a e?-wave resonance, which 
in turn tends to magnify the interference terms between 
p-d, and d-f partial waves, leading to a significant differ- 
ence between the thermalization times calculated with 
and without the inclusion of interference terms. 




FIG. 1: Averaged cross sections ({1112)} and (012) normal- 
ized to their value at zero temperature as a function of the 
final temperature of the mixture made of rubidium 87 in 
the two different internal states \5S 1 / 2 ,F = 2,m = 1) and 
\5Si/2,F = l,m = —1). The resonance shape results from a 
d-wave resonance. The quantity (0-12) is represented in dotted 
lines. The calculation of ((0-12)} made with interference terms 
(solid line) differs notably from the calculation where those 
terms are neglected (dashed line). 



IV. THERMALIZATION IN AN ISOTROPIC 
HARMONIC TRAP 

We have shown in the previous section that the relation 
between the collision rate and the thermalization rate is 
strongly affected by partial wave interferences. In this 
section we emphasize the role played by the confinement, 
and we derive the dynamics of the thermalization in all 
collisional regime [2(| . 

For the sake of simplicity, we consider the thermaliza- 
tion of an atomic mixture confined in a harmonic and 
isotropic trap: 



Ui{x,y,z) = -m % uj 2 {x 2 + y 2 + z 2 ). (15) 

By application of the method of averages, we derive the 
following set of six coupled equations starting from the 
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evolution of the square radius: 



d(r 2 ) l 
dt 

d(r ■ Vj)i 
dt 

dt 



2(r- v,. 



(r • Vilij)i, 



(16) 



The momentum and kinetic energy equations for the in- 
dividual species include a contribution from the inter- 
species collision term, expressing the fact that it is the 
total momentum and total kinetic energy of the system 
rather than the individual species' momentum and ki- 
netic energy that are summational invariants. These dis- 
sipative contributions to the species equations do not ap- 
pear in the equation for the radius since the collision is 
local and the number of atoms for each species is con- 
served. 

In the absence of inter-species interactions the set of 
equations (fTo) describes the monopole mode [13, [U, [HJ 
of both species independently. The collision terms involv- 
ing In do not contribute to the above equations because 
all the dynamic quantities are collision invariant for intra- 
species collisions. This is why there is no damping for the 
breathing mode of a single-species classical gas confined 
in a harmonic isotropic trap. In this particular geometry 
and in the presence of two species, the relaxation only 
occurs through the inter-species collisions whatever the 
collisional regime of a given species may be. 

The total number of inter-species collisions per unit of 
time Ti2 can be evaluated after the thermalization, once 
equilibrium has been reached: 



NiN 2 c(ai2) ( ixlo\IjJi (mi + m 2 )wiw 2 



r 2 72 V k B T 



f 



3/2 



(17) 



As the initial state of the mixture is assumed to be not 
too far from the final state, this expression gives a reli- 
able estimate of Ti2 during the thermalization process. 
The set of Eqs. ([TrJ)) is exact within the Boltzmann for- 
malism. To calculate quantitatively the relaxation driven 
by inter-species collisions, we make a Gaussian ansatz for 
the phase-space distribution function /j = /j(r, Vj,t) of 
each species, with the inclusion of a factor taking into 
account the space-velocity correlations: 



fi = Mi exp 



rrnvf + rriiOjfr 2 
2M 



(1 + rnmiT.Vi), (18) 



where Mi is the normalization factor. The effective tem- 
perature Ti and the parameter rji are the only time- 
dependent variables. Such an ansatz is inspired by the 
exact solution for the phase-space distribution function 
for the monopole mode, and provides a natural gener- 
alization of the local equilibrium distribution. For one 
species the Gaussian ansatz was shown to be accurate 
for investigating the damping of the coupled monopole- 
quadrupole oscillations in an anisotropic harmonic trap 
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FIG. 2: Evolution of the normalized moments Aj as a function 
of u>t in the collisionless regime ujtq — 20: Ai = Ai(i)/Ai(0) 
(solid line), A 2 = 2A 2 (t)/(Ai(0) + Aa(0)) (dashed line) and 
A3 = As(t)/ As(0). They account respectively for the differ- 
ence in the mean value of the square of the radii, the space- 
velocity correlations and the square of the velocity. 



from the collisionless to the hydrodynamic regime [17| . 
At the lowest order, the space- velocity correlations only 
enter the momentum equations (r-v,iy)j, while they can 
be neglected for the calculation of the term involving the 
kinetic energy (v?7y)i. The details of this calculation 
can be found in Appendix [Bl We find: 



(r ■ viii 2 )i 

(V?J12>1 



miu>l(r ■ vi)i - m 2 mj(r ■ v 2 ) 2 
m\Niu)iL02 f 

TOl(vf)l - 1712(^2)2 



mi N\ t 

where the time constants f and r are given by: 



1 
f 
1 

T 



[iUJ\L02 



.\0-12 



3 miwf + m 2 uJ2 (^12) 

8 H «<712» 



■r 12 



3 M (a 12 ) 



■r 12 . 



(19) 
(20) 

(21) 
(22) 



Combining Eq. flTBJ), and Eqs. (JTnj)-(E2), we obtain a 
closed set of six linear coupled equations. As a conse- 
quence, the relaxation in a trap does not correspond in 
general to a simple exponential decay. 

We demonstrate the physics of rethermalization using 
a specific example in which mi = m2 = m, uj\ = ll>2 = w, 
Ni = N2 = N and keeping a constant cross section. 
However, we emphasize that the conclusions we draw are 
generic. In order to follow the thermalization it is conve- 
nient to introduce the three quantities: 



Ai(t) 
A 2 (t) 
A 3 (*) 



vi)i(*)-(r.v 3 )2(t), 



'1/1 



(*) 



'2/2 



(t). 



We characterize the collisional regime through the 
quantity lutq with tq — Nr. The behavior we obtain is 
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FIG. 3: Evolution of the normalized moments Aj as a function 
of cut in the intermediate regime luto — 1, same notation as in 
Fig. ©. 




The time needed for thermalization is longer for a con- 
fined gas (see for comparison Eq. (|11J)). This is in agree- 
ment with the results of references]!, 0] ■ This factor of 
two accounts for the fact that the space of configuration is 
larger for a non homogeneous gas, and the thermalization 
affects both the space and velocity degrees of freedom. 

The decay exhibits oscillations when the collisional 
regime is such that lutq > 0.5, with an amplitude scal- 
ing as 1/wtq. In the intermediate regime (lotq ~ 1), 
the space-velocity correlations play a crucial role in the 
dynamics, and the decay of the Aj quantities is not ex- 
ponential at all as shown in Fig. [3l In the hydrodynamic 
regime (lutq <C 1), we find two separate timescales: a 
rapid one for the relaxation in velocity space, and a slow 
one for the relaxation of the spatial widths (see Fig. [4|. 
We recover in this limit the classical behavior of hydro- 
dynamics, with a short time needed to reach a local equi- 
librium and a longer time to reach the global equilibrium. 
The time needed to reach complete thermalization is pro- 
portional to 1/(cj 2 t ). In contrast to the collisionless 
regime, this time increases with the inter-species colli- 
sion rate. We have also performed numerical simulations 
based on molecular dynamics [23| to analyse the rether- 
malization in all collisional regimes and for arbitrary val- 
ues of atom numbers, atomic masses, and harmonic trap 
frequencies for the two species. We have obtained a good 
agreement with the predictions of our model based on 
the ansatz (IT%1) . 

In an anisotropic trap with different trap frequencies 
(Jf > ' y , one can derive a set of equations similar to Eq. 
(fT6|) . In this case the behavior of the relaxation is also 
influenced by the damping of the quadrupole modes [l7| . 



FIG. 4: Evolution of the normalized moments Aj as a function 
of cut in the hydrodynamic regime utq = 0.2, same notation 
as in Fig. 



reminiscent of the one of the monopole-quadrupole mode 
relaxation in a harmonic and anisotropic trap [17]. How- 
ever, the observables (r 2 ),;, (r • v<}< and (v 2 )^ and their 
equation of evolution are not the same. The different 
regimes depend on the dimensionless parameter ojt . 

For ujtq ^S> 1, the gas is in the collisionless regime (see 
Fig. [2]). We can recover formally the results obtained 
in the previous section, for a confinement of the gas in a 
box, by neglecting the space-velocity correlation term in 
the equation for A3, leading to an effective equation for 
thermalization that is identical to Eq. ([TT]) . In presence 
of the harmonic confinement, A3 obeys a third order dif- 
ferential equation according to the set of equations 1161 
In the collisionless regime, it can be recast, under our 
assumptions, in the form: 



dA 3 
dt 



A3 

T 



V. CONCLUSION 

In this article, we have investigated the thermaliza- 
tion of a non degenerate atomic mixture. In our analy- 
sis we have taken into account both the energy and an- 
gular dependence of the inter-species elastic differential 
cross section. We have derived the explicit form for the 
thermalization rate, which depends on the interference 
between the partial waves while the collision rate does 
not. In addition, we have shown how the dynamics of 
the thermalization is modified by the collisional regime 
in presence of a confinement. 
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APPENDIX A: COLLISION INTEGRALS 
CALCULATION IN THE HOMOGENEOUS CASE 

The global conservation of kinetic energy in a collision 
leads to E + = Si + £2 = 0. The symmetries of the 
collision integral permit one to rewrite the quantity E~ = 
Si — S2 in the useful form: 



liV / d 3 v d 3 v r d 2 n'v r ^- [v ■ (v P - v' r )] 



(AI) 

where we have introduced the quantity 2 
/iW,i)/iK,f) - /i(vi,i)/i(v 2 ,t). The Gaussian 
ansatz (|9]) for the distribution functions fa allows for an 

expansion of 2 2 around the final temperature Tf. We 
obtain 

x = --^Z- f^^MM / d 3 v d 3 v r d 2 n' 



k B T f 
da 1 



T r 



7 [v .(v r - v^fexp 



/ 



(A2) 



From the calculation of S~, one obtains the set of linear 
Eqs. (fTTj) with the explicit form Eq. (fl"2"|) for the relax- 
ation times r. The angular integration of Eq. (|A2|) leads 
to the quantity (Ti2(cx) defined in Eq. (|13p . From the 
expansion of the scattering amplitude in terms of partial 
waves ©, we can recast it in the form 



ai.i' sin Si sin Si> cos(<5; — Si>) 



0<2<i'<oo 

where we have introduced the dimensionless numeri- 
cal factors ai, v = (2 - 5j,|/)(2Z + 1)(2Z' + 1) f\(l - 
x)Pi(x)Pi' (x), (5/^' being the Kronecker delta symbol. 
The expression for the coefficients ai y explicitly con- 



tains terms describing interference between different par- 
tial waves. The properties of the coefficient a; are the 
following: a^i > 0, a;,z±i < 0, otiy — otherwise. 

Finally, we perform the Gaussian integration of Eq. 
(|A2|) . and we work out the expression for the relaxation 
time t of Eq. (fl"2f according to the notation of Eq. (jlljl . 



APPENDIX B: COLLISION INTEGRALS 
CALCULATION FOR A HARMONIC 
CONFINEMENT 

The calculation of the quantity (~v\li2) can be carried 
out in a way similar to the one used in Appendix lAl The 
terms in rji of the ansatz (|18[) do not contribute at the 
lowest order. 

To calculate the quantity (r • V1/12), we introduce the 
quantities A,; = rriiNi(r ■ ~ViIij)i and A ± = Ai ± A2. The 
global conservation of momentum ensures A + = 0. The 
first non-vanishing contribution to the linearized expan- 
sion of A~ is obtained by setting T\ = T2 = Tf in the 
ansatz (fTSl). One finds: 



(m - m)^^2 I d\ Q e- Mv o/ 2k - T f 

d 3 rd 3 v r d 2 rtv r ^-[r-(v r -V r )} 2 



exp 



2k B T 



f 



(Bl) 



The ansatz (1181) provides the explicit link between the 
space-velocity correlation moment and the parameters 
rji\ (r ■ Vi)i = 2>r]i(kBTf) 2 / (mnuf). This expression in 
combination with Eq. (|B1[) permits one to derive the 
explicit expression (|2ip for the relaxation time f of the 
space-velocity correlation equation. 
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